ANALYSIS OF LOSS OF COOLANT ACCIDENTS 


A Thesis Submitted 

in Partial Fulfilment of the Requirements 
for the Degree of 

MASTER OF TECHNOLOGY 


By 

RAKES H S HARM A 


to the 

NUCLEAR ENGINEERING AND TECHNOLOGY PROGRAMME 
INDIAN INSTITUTE OF TECHNOLOGY, KANPUR 

JANUARY, 1978 




This is to certify that this work on ’ihialysis of 
Loss of Coolant Accidents' has been carried out under our 
supervision and has not been ^submitted elsewhere for a 
degree. 




Dr. A. Sengupta 
Asst. Professor 
Nuclear Engineering and 
Technology Programme 
LIT Kanpur 



yn 

Dr. KT''"Srir^ 

Professor 

Nuclear Engineering and 
Technology Programme 
I IT Kanpur 


’■ ' ‘ * v'\ ^ T-f’ »-» 

„ ,..4 of . 

(he •a'vva^'-* ■ a.] 

; rTechsioiogY I 

::;,nce .athihe . - Y 

' ' - 'of 

' ;r:::.«o(Tcchn^ogy 



LIT. 


A 54»45. 

' 4i iMNi i|iM|k||H| 


■ 8 MAY 1P8 


NErP' 19Vs- |Yi _c.n A' 


ACKHOVniEDGMEHTS 

I am grateful to my guides Dr. K.Sii Ram and Dr. A. Sengupta for 
their guidance end cooperation throughout the course of this work. I would 
also like to record my gratitude to Mr.S.S. Pathak for his unstini-ted help 
whenever required. 

I am thankful to Dr. P. JIfegetnnaitecm. (bHED) for suggesting the 
problaa. 

I don’t want to thank Mr. K.Gargopadhyay and Mr. H. Kamick as 
it will reduce the importance and magnitude of their help. 

Mr. N.P. Jhaveri, Mr. P.S. Patidar, Mr. R. Jiyawan, Mr. A. Eajguru, 
Mr. S.C. Jethi, Mr. J.K. Deb and all other E-Block mates need; special 
mention for their respective help to my thesis work. 

I would also like to express my thanks to Shanma!, J. Singh, 

MM. Singh, R.I. Singh, E.S. Tripathi, R.M. Jha, Bharatiya, Ajodl^'a Prasad 
and mary others for their help, during the course of my thesis. 

Pinally, Mr. J.K. Misra deserves thanks for a magnificient job 
of typing and also Mr. Buddhi Earn Kandiyal and Mr. E.S. Tilpathi for u 
cyclostyling of the manuscript. 


Eakesh Shama 



G0MENT3 


Chapter 

I INTRODUCTION AND GUNSEAI RSUIEW OF THE 
PROBLEM 

1.1 Introduction 

1.2 Reactor Accidents and LOCA 

1.21 Reactor Transients 

1. 22 loss of Coolant Accidents 

1.3 Litemture Survey 

1.31 Theoretical Basis 

1.32 Experimental Methods 

1.4 Scope of the Present Work 

II THEORilTICAL ANALYSIS 

2.1 Introduction 

2. 2 Homogeneous Plow Model 

2.3 Separated Flow Models 

2.4 Determination of Quality 

III EXPERIMENTAL ANALYSIS 

3.1 Introduction 

3.2 Experimental Assembly 

3.3 Instrumentation 

3.4 Experimental Procedure 

IV RESULTS AND DISCUSSION 

4.1 Results and Conclusions 

4.2 Scope for Future Work 


Page 


1 

1 

2 

3 

4 

5 
7 

15 

16 
17 
20 
26 


31 

31 

36 

39 

40 
43 


REFERENCES 


45 



ABSTRACT 


The present work investigates the two-phase crittca.1 
flow rates in a loss of coolant accident of a PWR. An expe- 
rimental loop was set-up and flow rates, static pressures y 
stagnation pressures and temperatures were measured. A new 
type of quick acting va^lve was developed, made of burnable 
mylar sheets, which gave quick, clean breaks. A new ■ 

(indirect) technique to measure mass flow rate was developed 
by measuring changes in level. It was found that the 
technique needs refinement (by incorporating baffles, 
possibly) to cut out level changes due to vibration from 
the impact of the jet on the container. Unlike a,veraged 
measurement of flow this measures the flow continuously 
at all time instants after the break. 

An integrated approach to the study of equilibrium 
critical mass flows by the Moody and Pauske methods has 
been presented. Both Moody's and Pauskes* models were found 
inadequate for low quality flows. 
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tho cross sectional area of the break (m ) 

^ p 

the gravitational conversion factor (Kgm/kgf ) ( f t sec~ ) 

p 

the critical mass velocity (Kgm/ft -sec) 

the enthalpy of saturated water (Kcal/k.gm) 

the latent heat of vaporization (Kcal/Kgm) 

the enthalpy of the steam-water mixture (Ecal/Kgm) 

the slip ratio 

p 

the stagnation pressure (Eg/cm ) 

p 

the static pressure ('Kg/cm. ) 

the critical pressure (Eg/cm^) 

the entropy of saturated water (Ecal/Zgm °K) 

S.-S^ (Kcal/Kgm °K) 

I g 

the entropy of saturated steam (Kcal/Kgm °K) 
the entropy of steam f/ratcr mixture (Kcal/Kgm °K) 
the t hurst force (Kgf) 

tho specific volume of tho stoam water mixture (mr /Kg) 
tho steam quality 

the amount of the expelled water (Kgm). 



CHAPTER I 


INTRdPUCTIOR AND G-EIERAL REVIEW OP THE PROBLEM 


1.1 Introduction ; 

Energy crisis is one of the major problems the 
world is facing today. Many alternatives are being suggested 
to solve this crisis, like tapping Solar energy, more 
utilisation of hydli® and nuclear power, etc. As far as 
solar energy is concerned, the commercialisation has yet to 
take place, and most of the hydro power has already been 
tapped. Then comes, nuclear power, which is expected to 
effectively solve the energy problems for many years to 
come. To cope with the demand, arouiid 1000 nuclear power 
plants will have to be built by the year 2000. As the number 
of nuclear reactors will go up, there will be a strong need 
to study nuclear safety. The present work is an attempt 
to study an important aspect of nuclear reactor safety i.e. 
loss of coolant accident (LOCA) using an experimental loop. 

1.2 Reactor Accidents and lOGA ; 

In this section we will consider reactor accident 
in general and LOGA in particular. Reactor accidents can 
occur in many ways. This may be due to control rod failure, 
loss of flow due to pump failure, failure of reactor safety 
mechanism, etc. These accidents can be divided into two 
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major categories, firstly the design based accident (DBA) 
■which means that, kind of the event and severity are known 
prior to the accident, for example, loss of flow accident 
d'ue to p'ump failure etc. For such accident safety measures 
could be taken during the design of the nuclear power plant. 
The second category of accidents are all those events which 
cannot be predicted before the operation of the plant, e.g. 
refuelling accident, spent fuel handling accidents etc. 

Major consequence in all nuclear reactor accident is 
the release of large amovint of radioactivity to the atmos- 
phere due to the melting of the fuel in the reactor core. 

For a core melt-doim i.e, for maximum credible accident, 
either cooling system should fail or the heat imbalance 
should allow the fuel to reach a temperature of 5000°F. 

This would req'uire the reactor shut-down. 

A reactor safety study analysis^ has defined these 
two situations as 1. Loss of coolant Accidents and 
2. Reactor Transients. 

1.21 Reactor Transients? 

■■ ■' 

The term reactor transient refers to a number of 
events that require the reactor to be shut-do-^m. These 
range from normal shut-down for refuelling etc. to unplanned 
but expected events as loss of poi^rer to the plant from 
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utility transmission lines. The reactor is designed to cope 
with unplanned transients, but if all of them fail together, 
then the heat would be sufficient to eventually boil away 
all the cooling water and melt. the core. 

In addition to the above pathway to core melt, it is 
also possible to postulate the core melt resulting from the 
failure of reactor shut down systems following a transient 
event . In this case it would be possible that available 
cooling systems are unable to cope up and a core melt-down 
would result., 

1,22 loss of Coolant Accidents (LOCA) ; 

2 

The loss of coolant accidents are postulated to 
result from failure in the normal core cooling system. The 
water is at very high pressure in the reactor cooling systems 
and if a rupture occurs, in the pipes, pumps, valves or 
vessel, then a blow out would happen. In this case some of 
the water would flash to steam and blow out of the opening. 

The steam water mixture flows out of rupture until the 
pressure in the reactor and in the containment building 
become equal. This is known as blow down phase of a 
lOGA . For the first few seconds following the 

rupture, the fuel and cladding temperatures remain roughly 
constant, but with the onset of boiling about 5 sec. after 
the break, the temperature rises rapidly. For complex reasons, 
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having to do with the flow of ■water-steam mixture , out of 
the vessel, the cladding temperature, after reaching a. peak 
value of nearly 2000°F, drops down again and then undergoes 
a steady rise until EGGS (Emergency Gore Gooling System) 
intervenes. As soon as boiling occurs in the core, the 
reactor falls subcritical because of the negative void 
coefficient. 

For most accidents, i.e. core melt-down, lOCA must 
be followed by failures in emergency core cooling sys'’cems. 
But the principal exception to it 'is massive failure of the 
large pressure vessel. However, the accumulated experience 
with pressure vessels indicate that the probability of such 
accidents is very small and hence it does not pose much 
safety problems. 

1. 3 literature Survey ; 

Inspite of the great deal of work done on LOGA, the 
question of nuclear reactor safety has not yet received a 
satisfactory answer. One of the earliest reports in this 
direction was a study entitled, ’Theoretical possibilities 
and consequences of major accidents in large Nuclear Power 
Reactors carried out by Brookhaven National Laboratory in 
1957. Investigations into LOGA were also started around 


the same time. 
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Since the coolant in a majority of present day reactors 
is water, research in two phase flow has contributed a 
great deal both in terms of theory and experiments. .The theory 
has found application for no naqueo us coolants as well. The 
first full scale experimental facility lOFT (loss of fluid 
facility) was set up in Nov. 1964 by Idaho Nuclear Oorporation 
under the U.S.A.E.G. The other two set-ups were the ROSA-l 
in Japan and DHT-1 by ENBATOM at Ispra. 

1.31 Theoretical Basis : 

The final aim of all LOGA experiments is the prediction 
of the thermo hydraulic history of the corev This entails the 
measurement of Gq the inaximum critical flow rate at the rup- 
tured section, whicn will determine, along with initial condi- 
tions, the thermal hydraulic history in question. 

Several flow models have been proposed to determine, 

G-q, They can be broadly classified into two categories 
(i) Non equilibrium models (2) Equilibrium models, while 
equilibrium models have been widely studied, the non-equili- 
brium models have yet to establish themselves. Beloxir we give 
a brief account of non equilibrium models and a more detailed 
description of equilibrium models. 



The non- equilibrium theories assume. , the possibility 
of thermal non-equilibrium"' between phases due to finite 
evaporation rates. Henry and fauske treated non-equilibrium 
by imposing a condition of zero quality change between 
stagnation and throat. The vapour is assumed to expand isen- 
tropically and independently of the liquid fraction. The 
tendency of this limiting theory is "to over-estimate the nozzle 
flow rate. In general non-equilibrium fit better for short 
pipes, orifices and for stagnation qualities above 2% . 

The equilibrium theories, assume, a thermal equili- 
brium between the phases at the time of rupture. The first 
attempt in this diarection was made by Isbin et al*^ assuming 
a homogeneous flow model with the fluid properties averaged 
out. This model gave good results for large breaks, but for 
intermediate and small breaks, this model was inadequate. 

5 

fauske, H.K. , after a great amoinat of experimental work, 
assumed tha.t the two phase flow at the time of blow-down was 
annular with liquid in the annulus and vapour in the central 
portion of the tube. The model gave good results for all 
breaks but only for low pressures (0 - 360 psi). In his 
model, Fauske used continuity and momentum equations. 

C 

Moody proposed another mathematical model to calculate 
critical mass flow rate by using continuity and energy 
equations. A detailed account of these two models is given in 
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Chapter three. Many analytical models were proposed soon 

7 

after this work. MargoliSj, S.C. , and Eedfield devised an 
analytical model for blow-down experiments. This model was 
found inaccurate for describing experimental measurements 

O 

on the blow-down of simple tanks. They also studied reactor 
plant dynamics during lOGA. By assuming that the primary 
system was divided into three parts, rather than one lumped 
volume i.e. (i) pressuriser, (ii) primary coolant from the 
pump exit to the core inlet, he claimed that he got better 
results* 

1*52 Experimental Method s; 

The three full scale experimental facilities, as 
mentioned earlier, will be described first. Then we go 
over to other experimental works and application of computer 
codes in the study of lOCA. 

(a) LOFT, The loss of fluid integral test was designed by 

g 

Coplen and Tbarondo to study all the effects that were 
anticipated in a LOCA. The intention was to combine the 
possibility of studying nuclear, thermal hydraulic and 
structural proc esses occurring during a LOGA. Guret, : 
performed blow down tests to determine the influence of 
reactor vessel internal flow restrictions on the fluid 
decompression behaviour during a lOCA. Several different 
geometries were used to isolate the effects of hydraulic 



resistance and fluid inertia on the nature of i) core 
pressure drops occuring during "both the subcooled and 
saturated two phase blow- down iL)fluid flow separation , iii ) 
mass discharge from the system. Tests were performed with 
a hydraulically simulated scaled 10 FT core for breaks of 
various sizes in both the upper and lower plenums and it 
was found that the internal flow restrictions increase 
fluid phase separation and consequently increase residual 
water in the vessel after blow -down. The residual water 
was important as a source of steam for both core-cooling 
and metal-water reactions as well as for the calculation of 
core-f Idoding flow rates for EGGS. Using the same, Berta^^ 
studied subcooled pressure gradients in P¥R systems follow- 
ing loss of coolant. He used an approximation- method to dete 
mine the magnitude and duration of the pressure gradient. 

This method was an empirical formulation involving generali- 
sed system component dimensions and fluid properties. The 
results were compared with values obtained from the LOFT 
semi-scale blow-down test programme. The method consisted 
of evaluating a series of equations for each of the two 
plenums. The , solutions of the series of equations for the 
plenum closest to the pipe break gave the pressure rela- 
tive to the initial pressure in successive time steps. The 
time steps were computed by dividing the length of pipe 
between the break and vessel by the velocity of sound in the 



pipe. The same kind of information was provided for the 
other plenum in combination of two time steps, one being 
the previously defined one and the other calculated by 
dividing the core-length by the sound velocity in the core. 
Pressure gradients across the core are determined by 
differencing the pressure in the two plenum immediately 
above and below the core. 

Agreement between results obtained from the above 
calculation and lOPT data was good. The author claimed that 
the method was applicable t;;J a wide value range of parameters 
which include the dimensioiis of existing and proposed PWE's. 

12 

A.ET. Uahavandi et al. did a space dependent lOCA 
analysis for PWR's. The mathematical model used in the analy- 
sis was a very general and flexible one applicable to a 
large nximber of systems. The model allowed single, two 
phase and super-heat flow regimes & accommodated flow reversal 
between adjacent elements and choking flow at rupture points. 

The analysis employed a 3 dimensional approach (length, 
radius and time) in the reactor core, a two dimensional 
approach (length and time) in the reactor coolant loops and 
a lumped parameter representation in the secondary steam 
generating system. The fundamental equations of mass and 
energy for each element and momentum equation along axial 
directions were expressed and then these equations were 
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combined with l) energy balance equations and heat transfer 
correlations for the fuel rods and steam generator heat 
exchange processes, 2) subcooled, two phase and superheat 
state property correlations and 3) critical flow equations 
at the rupture elements. The resulting system of equations 
was solved numerically (SATAN Code) and the results were 
verified experimentally at loss of fluid test assembly- 
(LOFT). From this computer program the effects of hot leg and 
cold leg ruptures in pressurised water reactors, were deter- 
mined. Pressure, temperatures and flow rates were computed 
as functions of time at various locations in the reactor 
vessel with broken and unbroken loops. The lOCA was thus 
analysed in more detail. Calculations showed fair agreement 
with 10 FT tests, 
b) ROSA-I Assembly; 

The RDSA-I assembly was designed at Japan Atomic Energj?- 

Research Institute, its purpose was to study the thermo- 

hydraulic behaviour within the primary coolant system during 

13 

a lOOA. Makota Sabajima analysed lOCA with the help of the 
above assembly. In the ROSA-I assembly, the experiments were 
carried out by discharging water at high pressure and tempe- 
rature corresponding to reactor conditions from a pressure 
vessel, which was isolated from the heat up systems. Compari- 
son of the EOSA-I results with EBlAP-3, (a code used for the 
analysis of LOCA) revealed that results do not agree well under- 



certain break conditions, i.e, for hot leg rupture, the 
pressure and liquid mass calculated did not agree well 
with the experimental da^ta. 

The calculated pressure was lower than the data in 
the initial part of the transient and hi^er in the latter, 
finally it was determined that (i) the discharge coefficient 
used in the code as a para.meter could be correlated with 
the quality of the discharged fluid and (ii) the calculated 
liquid mass transient did not always agree with the experi- 
mental one when constant critical mass velocity was assumed. 
This discrepancy arose from the model dealing with the LOCA 
phenomena. 

c) DHT-1 loop; 

The DHT-1 loop (depressurization and heating transients) 
was built at the joint Research Centre Ispra, within the 
reactor safety program of EURATOM. The objective was to study 
transient heat transfer phenomena during depressurisation 
under P¥R conditions. The special loop DHT-1 were essentially 
one sub-channel and the upper and lower plenum of a P¥R. 

■ Friz and Riebold^^, with the help of this loop 
carried out a comparison between code calculations and blow- 
down experiments simulating a LOCA in P¥R. The test program 
was developed by G-AAA (Groupement Atomique Alsasienne 
Atlantique) at the same time but independently. This program 
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calculated the thermohydraulic history of the coolant within 
the core after a rupture in the primary circuit. The code 
calculated at different axial and radial positions within 
the core, the evolution of the coolant conditions and the 
temperature distribution within the rods and cladding of the 
fuel --elements. 

A series of experiments were carried out on DHT-1 
simulating rupture in the hot leg of a primary coolant circuit 
In the loop only three important components i.e. test section, 
upper plenum and lower plenum simulated by two additional 
volumeswere present. The loop did not contain pump or heat 
exchanger.- To ensure defined flow conditions prior to 
rupture, the whole loop was heated by trace heaters to induce 
natural connection upwanids in the test section. The measured 
parameters were pressure and temperature of the test section, 
whioh ara?« recorded on a multichannel recorder. 

The agreement between experimental and code calculations 
was found satisfactory, specially for large rupture cross- 
sections. It was- found that the experimental results agreed 
with calculations made on the basis of following assumptions. 

(a) Moody’s correlation with a contraction coefficient of 0.7. 

(b) Homogeneous flow model (i.e. slip ratio =1). 
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At the General Electric Company, a pipe rapture study 

15 

was done for the reactor primary circuit . It was c!esigr!,ed 
to study the mechanism causing failure and to provide 
improved capahilities for evaluating reactor piping systems 
to minimise the probability of failure. They made a) fati(p.e 
studies through bunch scale tests b) a failure probability 

x6 

study of actual reactor piping systems. Ishikawa,M= et al.'^ 

at JAERI have analysed blow-down behaviour of JPDR II 

(90 M¥(t) BWR), using the EELAP-3 code^"^. A series of tests 

were made on the number of nodes. y,'- ' . to 

IS 

obtain meaningful data. G. Gaspari et al. compared the 
scaled down experimental results with that from anadyrtic 
models and codes to predict the thermo -hydraulic history. 

The ez;^erimental measurements were carried out on a scaled 
down blow-down facility in terms of pressure, heated cha,nnel 
mass hold up and time to dryout under transient conditions 
determined by flow stoppage, break at different positions % 
power excursions. A blow-down code TILT was developed which 
has proved to be a useful tool for predicting the thermo- 
hydraulic history of the core. 

19 

Esposito studied the mathematical aspects of 
reactor blow-down and made a detailed analysis. To simulate 
lOCA, a large number of equations describing various thermo- 
hydraulic phenomena must be solved. Below a review is present 



of some of the existing computational methods used for this 
simulation. AppettjK.D. et.al. investigated the pressure 
oscillation that occurred during hlow-down in the pressure 
suppression system. It was carried out in the Marviken 
Nuclear Power Plant Sweden, from Aug, 1972 to May 1973. Ike 
object of the investigation was to determine the pressure 
oscillations occurring on the floor of the condensation 

chamber with ■ blow -down, also to determine the accompanying 

- 21 
mass flow composition of the mixture, ^3teinkamp; 'I * 

studied the analysis of the thermo hydraulic behaviour of 

P¥Ii plants in case of minor breaks in primary coolant 

systems. It was necessary to apply different approaches 

in this case, since the ruptured surface was so small that 

at no time during the course of the accident the core was 

completely exposed to the coolant. There was no refilling 

and rebonding phase as in case of large break lOCA. It was 

found that Pauske’s mode], holds good in this case. 

22 

Jenson . et.al. carried out blow-down heat transfer 
tests at Idaho Nuclear Corporation. The experimental data 
were presented in this report. The test was conducted at 
542 °P and 2260 psia (Reactor conditions), and a double ended 
cold leg rupture was investigated. The electrically heated 
core was used in the pressure vessel to simulate the effects 
,of a nuclear core. The flow resistance of the intact loop 
was based on core area scaling. Puring system depressurisation 
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core power was reduced from the initial level of 1,2 M¥, 
in such a manner to simulate the surface heat flux of lOPT 
nuclear fuel rods until., such time that departure from 
nuclear toiling occurs. Blow-down to pressure suppression 
was accomplished without EGG injection. The data was 
found useful. 

1.4 Scope of the Present Work ; 

The work done basically constitutes the first phase 
of an investigation into the loss of coolant accident, 

¥e have experimentally determined the two phase critical mass 
flow rate (Gc) of the coolant and compared it with the theo- 
retical values obtained using different flow models. 

Chapter two discusses., the two flow models used to predict 
maximum two phase flow rates, which form the theoretical 
basis of the problem. Chapter three contains an account of , 
the experimental work performed during the course of this 
study. In Chapter four we give the experimental results, 
and their interpretation. It also includes a few suggestions 
for further work. 



CHA.FTSR II 

THEORETICAL AMLYSIS 


2.1 Introduction s 

The general description of the loss of coolant 
accident and the need for analysis of such an accident in a 
nuclear power plant have been discussed in Chapter One, 

The present work simulates the loss of coolant accident in 
a PVR by causing a rupture in the pipe line. Due to the 
rupture in the pipe line the steam water mixture escapes and 
the behaviour of the escaping fluid is analyzed by consic eiring 
flow through an ideal nozzle. The purpose of this chapter 
is to study the two phase flow in this ideal nozzle. 

The maximum flow rate of a single component, two phase. 

mixture is of great practical importance because it restricts 

the capacity of the two phase flow rate. A number of models 

for obtaining the maximum flow rate of a two phase mixture 

23 

have been proposed. Among these are the workai of linning , 
Isbin , et al.^, Pauske ' ^ and Moody^ ■■. Of which 

the last two have been found suitable for blow-down experi- 
ments. The Linning.’ and Isbin models were based on empirical 
calculations. Moreover they assumed homogeneous flow model 
i.e. liquid and vapour are in equilibrium in the flow path 
(i.e. in the nozzle section) and the flow properties are 
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average of two phases. This approximation does not give 
good results for hlow-down experiment. Moody's and Fauske's 
models assume separated flow and they have been found sui- 
table under different conditions. The present work is based 
on- these two models. ¥e will first discuss briefly the 
homogeneous flow model in Section 2.2 and then go over to a 
detailed discussion of the separated flow models of Fauske 
and Moody. 

2,2 Ho mo ge ne o u s Flow Model i 

In homogeneous flow model it is assumed that a steam- 
water homogeneous mixture is flowing at a uniform velocity 
with the phases in equilibrium (Fo slip flow). The flow is 
steady and isentropic. The critical flow condition is defined 
under the restraints that the cha^nge in mass flow, rate with 
respect to pressure is zero at constant entropy i.e,, 




Maximum single phase, homogeneous flow rate - Homogeneous 
flow properties entering and leaving an ideal nozzle (entrance 
velocity is zero, total energy flow rates at entrance and exit 
are equal) are shown in Figure 2,1. 



18 



.ecus flow properties 


?ig. 2.1; Ideal. noszle and homogeneous li. p 
Continuity and energy equations are given hy, 


( 2 , 1 ) 


h 


0 


h + 


n 


£L 

U 



( 2 . 2 ) 


whe re , 

s = Mass flow rate/unit area, Ibm/seo-sq.f t 
h = Specific enthalpy, Stu/lbm 


U = Yelocity, f PS , 

Y = Volume, culic ft. 

subscript 0 refers to stagnation property 
Per isentroplc flew, h and Y are functions of P only 


h = h(s^, p) 

Y = Y(Sq, p) 


(2.3) 
( 2.4) 


where , 

Specific entropy KcalA^g' ^ 
o 

TO = Static pressure Kg /cm 
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Solving for G from (2,1) and (2.2) we have. 



2gc J ( 



) 


(2.5) 


Equations (2.3) > (2.4), (2.5) show that G i.s a function of 
P alone when h^ and are known. A maximum value of 6 
then must satisfy, 


dp ~ 

which leads to, 

" "®c f 

0 

Prom the . ' • thermodynamical consideratian, 

Tds = dh - 


( 2 . 6 ) 


(2.7) 


whe re , 

T = Temperature, degree P 

= J/144 = 5-403 ft^-lhf/in^-Btu 

if ds = 0 for isentropic process, it follows that 

Now (2.7) becomes, 

“ “ ^c ~ ^Iv^Sq 

which is the expression for critical mass flow of a single 
phase homogeneous fluid. 
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2*5 Separated Flow ModeLs ^ 

Both. I'auske and Moody developed expressions for the 
critical mass flow rate of a twq phase mixture assuming the 
vapour and liquid are separated in the nozzle with steam 
flowing in the outer region and the liruid flowing in the 
inner region. The other assumptions on which the models are 
based are; 

i) is en tropic flow, 

ii) the liquid and vapour are in equilibrium 
throughout the flow pathjf 

iii) average velocities of the two phases are different 
(i.e. slip flow is considered) 

iv) Plow is steady. 

Pauske'.s model is based on the continuity and momentum 
equations while Moody has developed a model which maximizes 
the theoretical two phase flow rate based on continuity 
and the energy equations. Different slip ratios are 
obtained in their models and appears that this discrepancy 
is due to the use of momentum equation in one model and the 
energy equation in the other, Note that, although the starting 
equations are different in the two models, the basic mathe- 
matical considerations to develop the critical mass flow 
equations are the same for both the cases. This point was 
unnoticed sofar and both the models will be described here 
together to understand this in a better fashion. 
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The continuity equations for vapour and liquid in a 
two phase flow are j 

G- X = a u ( 2.S) 

a a 

G-(l~x) = (l-oc) f f (2.9) 

where, 

G = mass flow rate per unit area of the homogeneous 

2 

mixture abbreviated 'mass flow rate' Ibm/sec-ft 
U = Average velocity ft /sec. 

'Z 

~ density Ibm/ft''^ 
a = Vapour volume fraction 
X = Quality of the mixture 


Subscript g refers to vapour and f refers to liquid. 


Kote that the left hand side of equation (2.8) and (2.9) 
describes the property of the mixture at the exit of the 
nozzle. 


In the absence of friction,, the momentum equation for 
such an isentropic two phase flow can be written as, 




Uf u^ 

(1-a) ^ — + a ^ 


+ 4 ^ = 0 

dz 


i.e. 4£ + £_ U (i/f ) = 0 

dz g. dz ^>m'^ 


( 2 . 10 ) 


= static pressure 


Vhe re , p 



22 


l/;p.^ = mean specific volume 


^ V 

a R 


4“ 


( I'^X,) ' 

Ti-ST 




( 2 . 11 ) 


Now the slip ratio K is defined aS; 
u 

K = ^ 

^f 

From (2.8) and .(2,9), 

K = 


( 2.12; 


X l-g V 9 

1-x a ^ 


( 2.13) 


Eliminating a between (2.11) and (2.13) 


1/f)^ (x, k, p) = -g^/[ 




dp 


], 


From Sqn. (2.10-) 


(1/p ) 

(x,K,P) = - g /[d --~5L] 


(2. .14) 
(2.15) 


'^Fauske dp "s 

Since Moody's model is based on the energy equation 
(and not the momentum equation), some additional properties 
of the flow have to be specified. These are v) entrance 
velocities are zero, vi) total energy flow rates at entranc 
and exit are equal. With these assumptions the energy 
equation takes the form, 


h 


- ^'^f 2g^ J 


) 


(2.16) 


where , 


h = specific enthalpy Btu/lbm 

Subscript o refers to property at stagnation state 
J = mechanical equivalent of heat^ 778 ft-lb^/Btu. 
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Properties at 
Entrance 

0 = 0 

h. 


B 


P 

Sg 

U, 


f - 0 
U_ = 0 



Properties 
at exit 

& 

h. ,h4^ 
o’ I g 

K 


P 

tJf 5 

\ 

X 


U 


Pig. 2.2: Ideal nozale and two-phase annular^, 
flow properties. 

Prom eqn. (2.13) we can write a as. 


a = 


1 


(2.17) 


where, V = 1/f = specific volume. 


Properties aalong the saturation line are functions of 
pressure only, 


Yf = 

Y^(p) 

®fg 

m 

II 



hf = 

hj(p) 

Y 

g 

= Y 
g 

(P) 

(2. 

18) 

Sf = 

s^(p) 


= h 
, g 

(P) 



11 

iifg( P ) 

®g 

= s 

g 

(P) 



s may be 

given 

in t erms 

of p 

and X 



OQ 

O 

II 

03 

(p ) = 

s^(p) + 

fg 

(p) 

(2. 

19) 


Finally equations (2. .8), (2. 9) and (2.16) through (2.19) 
may be combined to give the following relationship for G: 
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(x, K,p) = 




g 

J [h -h. -'^ 
c 0 -1 . s^g 

r [ 

i 

CQ 

1 

O 

CO 

2 



“-S -S_p 

1 0 f + _g 

O ^ .-r-kl 


’fg 




fg 


"fg 

( 2 , 20 ) 


where M refers to Moody’s model. 

If the stagnation entropy is known, then from 

-t-v, + p fbnth in ?auske's and Moody’s 
equation (2.19) it is obvious that G (botn ^ 

^ V function ol S and P only. Thiis & = G(K,P;. 

models) hecomes a function u 

• ■ -nnc,t o?atisfv the necessaiy 
Therefore, if G has a maxinun, ir must sat y 

conditions , 


(|S)p 


= 0 

= 0 


( 2 , 21 ) 
( 2 . 22 ) 


The condition (2.21) Sites the value of K at cntioal 
flow, from equation (2.20) and using (2.21) we get, 

1 . 






(2.23) 


+-0 nt critical flow for Moody's model, 
where = slip ratio at cntioal 

in Paushe's model, we have to use a condition which 

(0 01 ) to find out the slip ratio 
is equivalent to equation (2.21) to lin 


namely > 
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dE 


= 0 


. (2.24) 


Eauske, in his model, assumed that at critical flow 
the pressure gradient attains a finite mazimum mlue. 
Now, differentiating (2.10) with respect to K 


-A 

dK 


G d& d 


+ r- ^ -rr (!/?„) + ^ (l/f„) = 0 


g dK dz 
°c 


I m 


g dK dz 
“c 


Using Kauske * s assumption and equation (2,21) we get 


dz dK ° 

i.e. ^ (l/fm) = Const. 

Fauske assumed this constant to he zero and thus employed 
equation (2.24) to calculate the slip ratio for critical 
flow. 

Thus .from (2.24) and (2.14), we get, 



where Kj, = slip ratio at critical flow for Fauske 's model. 

KTow, using eqn. (2.25) iu (2.15) and eqn. (2.25) iu 
eqn. (2.20), G becomes a function of p only. To get the 
critical flow rate, we now use the condition (2.22). For 
this a numerical technique is employed and the prescription 


is as follows ; 
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Assume a value of quality x, and then calculate & for 
various p, the point of inflation is noted. The quality is 
changed, and the calculation of G- is repeated for- various 
p to get another point of inflation. Thus the maximum flow 
rate is obtained. 

¥e will describe the indirect experimental proce- 
dure to measure quality of the two phase mixture. 

2,4 Determination of Quality s 

As we have seen that earlier methods can predict 
the critical flow rate for steam water mixtures from known 
static pressure, enthalpy and quality conditions, but 
these methods cannot be applied to the vessel blow-down 
calculations because the quality of steam water mixture 
cannot be experimentally measured. 

The present work is based on an approach ducj. 

24 

ta Shiba and Guret to determine the qualit 3 r of 

the mixture using pressure data available from experimental 
and two phase flow models. With the quality information, 
the critical mass velocity can be calculated and integrated 
through the process of the blow-down and multiplied by the 
cross sectional area of the flow path to determine the 
amount of expelled water. The accuracy of the critical 
mass velocity calculations can be determined from the 
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comparison of the calculated and measured values of 
expelled water, during blow-down tests. 

Thru-st and Gritioal Flow Relatio nship '; 

Prom the earlier developments, we see that when 
appropriate value of K is substituted in equation for G, 

G becomes a function of p and x in both Pauske's and 
Moody ’ s models , 

G = G(p, x) (2.26) 

Prom the simulated experimental assembly we get the stagna- 
tion pressure within the vessel, during blow down. The 
experimental pressure and thrust data in conjunction with 
analytical techniques can be used to determine the quality 
required to predict the critical mass velocity for vessel 
blow down. The thrust applied to the cross-section of flow 
path is 5 

T = P X A (-2.27) 

where P is the stagnation pressure, and A is the cross- 
sectional area of the break. 

Now,, from, a simple force balance, the thrust is the 
summation of the momentum expulsion rate and the pressure 
force, 



1 

A 


( 2.28) 
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where p is the static pressure, and ^ (the mean density) 
is , 


fm 


[x Vg + (1-x) K 


[x + 



( 2. 14 ) 


From equation (2.28), (2.14) and equation (2.27), one can 
easily obtain G- as a function of x and p. Uow, from this 
equation of G- and equation (2.26), (Moodj^'s Pauske's G), 

X can he eliminated and G becomes a function of p only. Then 
the application of the condition (2.22) becomes straight- 
forward. But the elimination of x my become cumbersome 
and hence we tried to get rid of this by suitabljr using 
experimental results in eq. (2.2S). Prom the experimental 
data, amount of water expelled is known for a given period 
of time if the differential pressure between the system 

pressure and outside presstire is assumed to be a constant. 

t 

¥ = A ( G- dt = A G t (2.29) 

b 

The critical mass flow is obtained from the corresponding 
maximum amount of water expelled. 


¥ ~ 
max 

Then equation 


A 


4 t 

(2.22) for critical. 



flow becomes, 


( 2 . 50 ) 


( 2.51) 



where , 

= static pressure corresponding to critical 
flow (experimentally measured). 

Gq =r critical mass flow rate 

Ifote that, ^ ^ is a function of pressure. Putting 
the mlue of G from (2.30), in (2.31) <? is calculated. 

^ * D3. 

Then from (2.14), the quality x is known. Once the quality 
is thus experimentally determined, it can be substituted 
either in Moody's or Fauske's equation (2.26). Then by 
virtue of (2,22) the theoretica.1 G can be calculated. This 
theoretical value can be compared with experimental G^. But 
this procedure has a big diawback in the assumption that 
the differential pressure between system pressure and 
outside pressure will remain constant. 

21 ? 

This driving differential pressure is very important"^ 
because it is used to overcome the acceleration, frictional 
and other pressure losses of the escaping fluid up to the 
point of break. In a system, where the volumes inside and 
outside are very large compared to the size of break, the 
differential pressure can, be assumed constant. But the 
experimental assembly, which Curet has used, the pressure 
drops to its 50% value with in one second, this assumption 
is not valid. So the critical mass velocity, which we will 
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get by this model, -will just be an average escape velocity. 
We have tried to modify this system, by assuming the 
escape velocity constant for a fraction of second and then 
tried to use the same analysis. Moreover, we are 
indirectly measuring the critical mass velocity continu- 
ously and recording it on the storage oscilloscope-. We 
feel that this assumption will hold good for calculating 
the critical mass velocity in our approach. 
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EX PERT MBFTAL ANALYSIS 


3*1 Introduction ; 

Many experimental assemblies have been built to study 
the loss of coolant accident. The present experimental 
work has been performed to check the critical flow rate models 
i.e. choked flow models. The prediction of critical mass 
flow rate is important for analysing reactor safety because 
this -determines the maximum escape rate of a high pressure 
fluid from the system in which it is contained, which in 
turn fixes the pressure, that the containment has 
to with stand. 


3 . 2 Experimental Assembly ; 


The schematic diagram of the assembly is shown in 
figure^l) . It consists of a test section, with upper plenum, 
lower plenum, a water pump^ tanks to collect the fluid, quick 
acting valves and a flowmeter. Each of these will now be 
described in detail. 

(a) T est Section ; 


The test section consists of a mild steel pipe, of 1 in, 


internal diameter iVs iiich. outer diameter. It has five 


LIT. m 

cmmx-L iimMi 

Ace. Sfe. 



TEST SECTION 















P°PP"^ i rnnstoot 



Transient Water Level Measurement 




openings for fixing thermocouples, made of copper-constantan. 
The test section is heated externally hy a heater ■which was 
fabricated in our laboratory. Originally, we had planned 
to simulate the heat flux of a 220 MV/e P¥R reactor in our 
set-up, but due to experimental limitations, this could 
not be achieved. A heater of 3.2 Kilowatt was fabricated 
with a 20 S¥G twisted (double) nichrome wire. Initially/, a 
double layer heater was tried, we tried to insulate the 
layers by glass -fibre, and then by thick mica sheets, but 
due to frequent short circuits, this idea was given up. The 
nichrome wire -was wrapped on the test section. The major 
difficulty in wrapping of the resistance wire was the 
separation between consecutive wires, as less than adequate 
separation resulted in a short circuit. 3y trial and error 
method we found that the separation between two wires should 
not be less than 3 i’mi. Also, a -mixture of white cement and 
asbestos powder was filled in the space between the wires. 
Subsequently, thermal insulation 'wa.s provided by pasting 
asbestos and white cement mixture around the heater. The 
higher heat flux can be achieved by the use of an induction 
heater. 

( ) Unper and lower Plenum ; 

The upper and lower plenums of the reactor were 
simulated by using two 3'* inner diameter mild steel pipes, 
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which had openings for q.uick acting valves and couplings 
for pressure transducers,.. The plenums were fahricated in 
our central workshop.. The, holes were drilled on the plenum 
surface, where counter part of quick acting valves were 
welded.. The water- expelled was collected in aluminium tanks 
of dimensions 1' x 1-p x l|’- *' . These tanks were fahricated 
in the central workshop. The lower plenum has two tanks, 
one for double ended break and another for single ended 
rupture. The upper one has only one tank for single ended 
rupture.. 

(c) Water Supply Syst em; 

The water was stored in a tank 6' x 4' x 4' from which 
a four stage, high pressure, centrifugal pump forced the water 
to the test section. Specifications of the pump and motor are 
given below. 

Pump ; BE - 3141 Model A 0214GK 

Capacity 3 gpm at 300 psi (ma-ximum) 

Mo to r Crompton Parkinston, Bombay 
1450 PPM, 4 Phase 440 Y 
7.5 A, Mesh, 5 H.P. 



(d) F low-meter s 


A rotameter was attached to the pump in order to 
measure the flow rate of the water. It was a variable area 
flow meter with a capacity of 10 gpm. 

(e) Pressure P robe s; 

The pressure is measured at four places in the loop. 
Bourden gauges were placed at lower and upper plenums respec- 
tively to measure the systems pressure. The pressure 
.transients were measured at two places, one inside the lower 
plenum and another just at the throat of tho break. These 
were measured by two strain gauge type pressure transducers, 
manufactured by Consolidated Electrodynamics, Pa^dena, 
California (a subsidy of Bell and Howell). The specifications 
were as followss 

Humber Eange Max. Temp. Hat.Freq. 

CEC 4-313-0001 0-300 psia 250°F 15 SHz 

4-312-0001 C-I50 psia 250°F 15 SHz 

The transducers did not have connectors and cables. The 
connectors were fabricated on a printed circuit, in the 
printed circuit lab. and precision shop. 

T) Quick Acting Valves ; 

One of the important work done during this thesis was 
to make quick acting valve. This valve was needed to 
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simulate the rupture in a pipe. In general the rupture time 
in a nuclear power plant is less than a millisecond and it 
is difficult to have such a fast valve. In addition to 
these aspects 5 the valve must operate at high temperature 
( ~ 100°C) and high pressure ('•^ 200 psi). InitiallTr, we 
started x-rlth a spring device, hut it had imny limitations, 
like the high time constant (2^50 ms). Moreover, it was a 
difficult task to make it leak proof. Ouret and Shiha^^ 
used two discs, which were ruptured by pressurizing the air 
between them. This device has man 3 r disadvantages. For exam- 
ple, the pressure which is induced externally gives a shock 
wave which travels opposite to the internal pressure shock 
wave. Secondly a big compressor will be required which will 
create 'many more problems. Accordingly we tried few 
different approaches for creating the rupture. In one, the 
rupture disc, was burnt by an electric arc. A high voltage 
was applied between rupture disc (aluminium foil) and an 
electrode, by an induction coil (15 KV), But it did not work, 
as the aluminium foil took a long time to burn and secondly 
i-t could not sustain high pressure. Another attempt was 
made to construct a device which could be controlled by a 
magnetic field. It was found, however, that a very high 
magnetic field would be needed to operate tnis device and 
the idea was accordingly given up. 



36 


Finally, we tried plastic discs, which could sustain 
high pressures and temperature. Ordinary plastic sheet was 
used, hut at high temperature it melted down and the attempts 
to insulate them thermally went in vain. This led us to use 
mylar sheets (a type of plastic which can sustain high tempe- 
rature ('---'150*^0) and pressure (^■^200 psi)). Two mylar sheets 
with a resistance wire in the form of a coil of a fixed 
diameter sandwiched in between formed the rupture disc. The 
resistance wire was heated externally and it took less than 
Ip sec. to rupture the disc. Thus a new technique was developed 
to simulate rupture in a pipe i 

5 • 3 Instrumentation ; 

The instrumentation was designed to measure the follow- 
ing variables. 

(a) Stagnation pressure 

(b) Pressure transients 

(c) Temperatures, and 

(d) Flow rates 

( a ) Stagnation Pressure ; 

The stagnation pressure i.e. the pressure inside t.ie 
loop was measured at two places, (one near lower plenum, 
another near upper plenum) by standard Bourden gauges, which wer 
available in our laboratory- 



("b) Pr e ssure Trans ienta ; 


The pressure transients were measured for both static ' 
and stagnation pressure, as soon as the rupture occured. 

These transients were measured by two pressure transducers. 

The static pressure was measured just at the throat of the 
break and stagnation pressure transients were measured at 
the lower plenum.. The static pressure was recorded on a 
storage oscilloscope (BGIl, 768 3). The stagnation pressure 
transients were recorded on a multi-channel recorder 
(Encardio-rite, 4 channel strip chart recorder). The oscillo- 
scope can not take differential input, so one differential 
amplifier was designed"^, for the pressure signal. The 
manuals for these transducers were supplied by Bell and 
Howell. 

( c ) Temperatures 

There were sewen locations- to measure the fluid tempe- 
rature. Pive places were in test section and one each for 
lower and upper plenum. Copper constanian thermocouples were 
used. The temperature was recorded by three recorders at 
three locations. Other locations were set idle as we did 
not have multichannel temperature recorder. At one location, 
i.e. at lower plenum, the temperature was measured by a 
milli--volt potentiometer, which could read variations of 
0.1°C. At the second location i.e, middle of. test section 



it was recorded by a single channel temp, recorder (Leeds and 
Forthrup 5 1 channel Time constant 2 minutes). At third location, 
i.e, test section, it was recorded on encardio-rite multi- 
channel recorder. 

( d ) Flow-rates 

Two flow-rates were measured. First one, is the 
pumping flow rate i.e. the amount of water flowing through 
the loop. This was measured by variable area flow meter. 

The other more crucial flow rate is the, critical mass 

flow rate. One indirect method, which was used by Shiba 
24 

and Guret , is to measure the amount of water collected in 
certain amount of time, but it has one drawback, i.e. they 
assumed flow rate to be constant throughout the process of 
flow down. In .fact flow rate does change and we developed a 
method to measure flow-rate continuously as the pressure 
varies. An attempt was made by using strain gauge, but it 
did not work, as the vibrations due to pressure and burst were 
too much and we got spurious signals. The other technique is 
to use load cell, which can record change in amount of water, 
but the water tank cannot be mounted on the top base in.dia) 
of the load cell. The technique employed in our work was 
simple but very effective. The principle is to measure resis- 
tance of water as it changes. This was done by putting two 
copper wires inside a glass tube, which is attached to the 
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water tank as shown in Figure 3*3 • A constant current (d.c.) 

supply was given to the copper wires. The voltages across 
the wires changes with the -water level. This change in 
voltage was recorded on the storage oscilloscope (ECU, 768 S) 
The output was of the order of 0-5 volt/mm change in water 
level (current 20 mA). The constant current supply was 
designed and fabricated in our laboratory^ 

3 • 4 Experimental Procedure ; 

Initially, cold water was filled in the thermally 
insulated loop. The water is circulated in the loop in 
equilibrium with the help of 4 stage pump. The heater was 
switched on and the water was circulated, till the tempe- 
rature went upto 100°C. The oscilloscope, Encardio-rite 
recorder temperature recorder and the various supplies were 
switched on. The temperature at lower plenum was measured 
with the help of milli-ohm-potentiometer continuously. As 
the temperature reached 100°C the mylar discs were burnt 
with the help of a high current ( 50 Amp). The oscilloscope 
was kept in a single sweep mode and the signals triggered 
it automatically. The static pressure transients and water 
level change, were recorded simultaneously on storage 
oscilloscope i The temperature signal near the lower plenum 
and stagnation pressure transient were recorded en encardio- 
rite strip chart recorder. 
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EESULTS AND DISCUSSION 

4.1 Results and Conclusions ; 

As mentioned in Chapter III, the values of the stagna- 
tion and. static pressures, and water level were measured for 
fo\ir different pressures. The temperature transient was 
measured only for 100 psi. Since no significant transients 
were observed the temperature measurement was discontinued for 
other pressures. The pressure transients (static and stagna- 
tion) are shown in Fig. 4.1 and 4.2. The lack of temperature 
transients arise mainly due to the low heat fluxes involved. 
Table 1 summarises the readings obtained together with the 
values calculated from Moody's and. Fauske's models. 


Table 1: Calculated and Measured Values of 0^ 


Set 

No. 

Stagnation 
Pressure 
Kg/ cm^ 

Static 

pressure 

Kg/cm^ 

Measured 

&c 

Kg/m^ec 

Calculated 

Moody's Fauske's 

Model Model 

Kg/m^“Sec Kg/m^-sec 

1 

3.52 

1.C7 

14702 

8480 

8697 

2 

4.40 

1.09 

14260 

3885 

3035 

3 

6.10 

1.10 

13240 

3388 

1614 

4 

7.80 

1.27 

12220 

2661 

1271 


The first 

calculation 

(Moody's) 

requires 

stagnation 


properties for the steam— water mixture and also the specific 
enthalpies for the pure gas and li(3.uid at the static pressure. 
Since, in o\rr case the inlet mixture contained no steam, the 
stag^nation enthalpy was that of the liquid alone. The property 




Fig.4-1 Stagnati 




Time ms 

Fig.4-2 static Pressure Transients 
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values were taken from the ISI stccmi talles^^. The calculation 
ot hy iJ.a,UEke s mod cl was socicwhat more tedious since various 
differentioas with respect t’^ pressure had to be calculated. 
This was done by the standard method of dividing ordinate 
differences by abciwssa differences (iis^/AP). The values, 
agruin, being obtained from the steam tables. 

From the results above it seems that both Moody's and 
Fauske’s models predict -values which are significantly 
different from the experimental ones.. It is difficult to 
say in the present case which, if any, of the two models 
apply- 

In our view the following are probably the major 
factors which led to the above discrepancies and which need 
to r'.-C i;ifii-d in future experimental work on the present 
s u t-up 

1. The main source of error G.:.Gms to lie in the measured 

V'llui r> or flow „ rate. The rather short time spans 
involv..-d 50 ms) made it necessary to measure the level 
ch'’!.nfrw r'‘th’-r tha.n the averaged flow rate th8.t has usua.lly 
b I,, n nu n.'d by other workers in the field. The meaoure>-. 

nun-it vni;- affi.cted by the impact of the ejected fluid end the 
r corded 1 vc.l change was significantly more than the actual 
ch';,ngw in l.vca. This was also tho reason why pressures 
hig’ru-r than 100 psi were not considered. An accurate 
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exp- rinent.nl arrangement would probably require baffles to 
break up the impacting fluid, ideally a very fast (open-shut) 
quick acting valve with variable ’open’ times would solve the 
problon. 

Fauske's model uses the experimentally measured flow 
rates to calculate quality and the discrepancy between the 
two c".lculatcd values at the higher pressures (6.1 and 7.3 Kg/cn^ ) 
can probably be traced to high values of the measured flow rate. 

2. At the inlet (stagnation point) the two phase mixture 
contained only liquid since a boiler was not available. This 
led to vury low qualities after the break (static point). 

3. The lov.' heat fluxes obtained from the heater prevented 
the systCim from reaching high temperatures. This again con- 
tributed to low qupility. 

4. 7h.. total water in the system was quite small ( 7-B IttftC.) 
Thus the period for which critical flow actually existed was 
ouitG short, of the order of 10-20 ms. This was a big handica.p 
since the fastest recorder avo-iloble had .a minimum tine base 
of 10 mo. Both static pressure and temperature me.asurements 
cuff rud fro'-^ this drawback. 

5. Th.-,. short length of the break meant that the assump- 
tion of Liuailibriur .1 between phases was seriously in error, 
both t:i- Fauck-- and Moody calculation have been vitiated 

by this aliortcoming. 
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6. Fauske’s and Moodjr's models are not applicable to low 

quality flows. As mentioned earlier the output quality -was 
low, the maximum being ?.S7 and it is doubtful whether the 
models actually apply, however in the absence of any other 
these two had to be used for purposes of comparison. 

4 , 2 Scope for Future Work : 

(a ) Sxporimontal ; 

One of the major problems in all LOCA simulations has 
been the extreme combination of conditions involved. For 
•„'XamplC‘ the high licat flux in nuclear reactors is quite 
difficult t ) get in a laboratory environment. Besides, 
critical flows in simulations exist for a very short time, 
of thu order of milliseconds whicb makes sophisticated 
inatrumentrtion necessary. 

No existing experimental assembly actually simulates 
n LOCA in ' real nuclear reactor. So a nuclear reactor meant 
specifically for LOCA exp.'riments is being constructed in 
the United States. 

’'/hllw an attempt has boon made in the present sot-up 
'.chi .‘VC- reactor conditions, several improvements ne>.;d to 
b',, "I'ldi, . T.h.osiO ca.n bo summarised as follows: 

( ) Vo »'!atch h'.at fluxes an induction heater can replacva 

th r(.i;iDt'‘nc;. wire he-ator being used at present. 
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(ii) A tio.-.t exchanger enn ho xirovided in the loop to 
siniulatu th_ primry heat exchanger in the reactor. 

(iii) DinenaionlcsB parameters like Reynolds and Prandtl 
numb.^ra can made eoual to those, in reactors . 

(iv) A aurg’. tank can be provided to alter systems pressure 
without cha.nging flow rates. 

(v) A multichannel recorder which can record fast 
tr'insicnto C'-n be attached to the system. Better still, 
the signals ca.n b ; recorded on -a real time data acquisition 
ay at '.I likv, the IBM 1800. This will give a detailed picture 
of the t r .nnients. 

( b ) flu. ,)re tlcal : 

Th . two :n?ii r I)rosont equilibrium models use either 
tii. rgy (M<a;/’y) or the moaontui (Pausko) equation along with 
thw c ntinuity equation to riaxinisc the two phase flow rate. 

A gi.nwrii aipro ich using the continuity, momentum and 
nu rgy '-auai.i n with the three independent variables, 
pi*<..: 30 ur , r:li] rritic and quality would probably give a 
b'.’t t .r 'i )'‘’.cl, 

Th'. equilibrium models, however, are clearly inadequate. 

In n n no c.an the two phases bo said to be in equilibrium 

29 

t'l.. brink. This means th-t non-equilibrium models 
hav- t'* b« '’'ove loped and experimental results compared with 
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